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Abstract 

We present the results of a new evaluation of the running fine structure constant 
a at the scale of the Z mass in which the role of the e + e~ annihilation input data 
needed in this evaluation is minimized. This is achieved by reducing the weight 
function M|/(s(M§ — s)) in the dispersion integral over the e + e~ annihilation 
data by subtracting a polynomial function from the weight function which mim- 
ics its energy dependence in given energy intervals. In order to compensate for this 
subtraction the same polynomial weight integral is added again but is now eval- 
uated on a circular contour in the complex plane using QCD and global duality. 
For the hadronic contribution to the shift in the fine structure constant we obtain 

(5) A 

Ao^ ad = (277.6 ± 4.1) • 10 . Adding in the leptonic and top contributions our final 
result is a(M z y x = 128.925 ± 0.056. 
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1 Introduction 



Currently there is a great deal of interest in the accurate determination of the running 
fine structure constant a at the scale of the Z mass |T|, |3|, |J. The value of a[Mz) is 
of paramount importance for all precision tests of the Standard Model. Furthermore, an 
accurate knowledge of a(Mz) is instrumental in narrowing down the mass window for the 
last missing particle of the Standard Model, the Higgs particle. 

The main source of uncertainty in the determination of a(Mz) is the hadronic contri- 
bution to e + e~ annihilations needed for this evaluation. The necessary dispersion integral 
that enters this calculation has in the past been evaluated by using experimental e + e~ an- 
nihilation data. Disparities in the experimental data between different experiments suggest 
large systematic uncertainties in each of the experiments. In order to reduce the influence 
of the systematic uncertainties on the determination of a[Mz) the authors of []3], |^| have 
added theoretical input to the evaluation of the hadronic contribution to a{Mz)- 

Global duality states that QCD can be used in weighted integrals over a spectral func- 
tion if the spectral function is multiplied by polynomials (i.e. moments) but not if multiplied 
by a singular function such as the weight function H(s) in the present case (see (fj)). Never- 
theless local duality is expected to hold for very large values of s, i.e. Im II(s) ~ Im n^ CD (s). 
The authors of || use QCD perturbation theory in the form of local duality in the region 
above s = (1.8 GeV) 2 for the light flavours. The authors of U use perturbative results for 
energy regions outside the charm and bottom threshold regions. In the respective threshold 
regions they use renormalized data where the renormalization of the threshold data from 
each experiment is carried out by comparing with QCD perturbation theory outside of the 
threshold region. The authors of ||, [|| differ in their treatment of the data. Both of the 
evaluations suffer from assumptions on the nature of the systematic uncertainties of the 
data. 

Our approach is quite different. We attempt to minimize the influence of data in 
the dispersion integral over the whole energy region including the threshold regions. The 
essence of our method is to diminish the size of the weight function in the dispersion 
integral by subtracting a polynomial weight function which mimics the weight function in 
given energy intervals. In order to compensate for this subtraction the same polynomial 
function is added again, but now its contribution is evaluated on a large circular contour 
in the complex plane where perturbative QCD can be safely employed. 

2 The method 

As is well known (see e.g. jl]) the hadronic contributions to the effective fine structure 
constant can be expressed in terms of a weighted dispersion integral over the total e + e~ 
hadronic annihilation cross section R(s) or, equivalently, over the imaginary part of the cor- 
relator of two electromagnetic currents. Our normalization is such that R(s) = 127rlmil(s) 
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where the four-transverse piece of the current- current correlator H(s) is defined by 



i J (0\jr(x)f^(0)\0)e iqx d 4 x = {-g aP q 2 + q a qpMq 2 ) 



(1) 



The hadronic contribution to the fine structure constant at the scale M z is determined by 
the dispersion integral 



where the weight function H(s) is given by 



The physical threshold for the light (u, d) -quark currents lies at So = 4m^. The physical 
threshold for the strange quark is nominally higher but can be lumped together with the 
(u, (i)-quark threshold for our purposes. Henceforth we shall therefore refer to a common 
light quark threshold of So = 4m^. for all three (u, d, s)-quarks. For the heavy quark 
currents we take the masses of the lowest vector quarkonium states as threshold values, 
i.e. the relevant thresholds values are given by m% and m\. 

The usual procedure to evaluate Eq. (0) is to substitute the experimental cross section 
into (0) up to some high momentum transfer value si, and from then on to replace R(s) 
by 127rImn^ CD (s) hoping that QCD furnishes an adequate description of the data above 
this momentum transfer value. Since the weight function in the dispersion integral far 
from the Z-pole is essentially given by 1/s, the phenomenologically determined low energy 
part of the dispersion integral dominates over the perturbatively evaluated high energy 
part making latter substitution quite safe. We want to make two remarks prompted by 
this observation. First, one tests very little of perturbative QCD in such an evaluation 
even if the data were perfect. Second, a serious drawback of such an evaluation is that 
the total hadronic e + e~ annihilation cross section is beset with large systematic errors 
that will directly feed down to the evaluation of aiMz) and make such a calculation quite 
unreliable. One may thus hope that the inclusion of additional theoretical input on the 
strong interactions may reduce the resulting error in the integral of Eq. (|2|) . 

In this paper we shall extend and elaborate on a technique proposed in || that allows 
one to substantially enhance the contributions of perturbative QCD to the evaluation of the 
dispersion integral in Eq. ([|). Let us explain the method for the case of the light u-, d- and 
s-quarks. We split the region of integration into two parts, one from the threshold so = 4m^. 
to some large value s\ where perturbative QCD is valid, and the other from s\ to s — oo. In 
practise we shall further subdivide the interval from s\ to s = oo into smaller intervals but 
this need not concern us here. As the weight function H(s) is not an analytic function QCD 
cannot be directly employed in the region s < si, the low energy part of the integral in 
Eq. (0). The concept of our approach consists in constructing a polynomial function Pjy(s) 
which approximates the function H(s) in the energy interval Sq < s < Si. The polynomial 




(2) 
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will be determined by the method of the least squares with possible additional constraints 
which will be discussed later on. One then adds and subtracts the polynomial function 
from the weight function. In the subtracted piece one achieves a substantial reduction 
of the influence of the data and their errors on the evaluation of the dispersion integral, 
whereas the added piece can be evaluated by using perturbative QCD. 

Accordingly we now proceed to write down an identity for the low energy part of the 
integration in Eq. @ by adding and subtracting the polynomial function Pjv(s) in the 
integrand. One obtains 

[ S1 -lmU(s)H(s)ds = r -lmU{s)(H(s)-P N (s))ds+ H -lmU(s)P N {s)ds. (4) 

Since the weight function Pn{s) in the second integral on the right hand side of Eq. (|j) 
is analytic one can replace the integration on the real axis by an integration on a circle in 
the complex plane such that one has 

j Sl -lmU(s)P N (s)ds = — — / U(s)P N (s)ds. (5) 

Js 7T 2lli J\s\=si 

Upon substituting (|) into (£D one obtains an exact sum rule if full QCD were used to 
evaluate the different contributions. However, since the full QCD expression for the current- 
current correlator is not known, we have to use the sum rule in an approximate sense. As 
explained before the first integral on the right hand side of Eq. @ will be evaluated using 
the data set given by 0. The small weight factor (H(s) — P/v(s)) will significantly reduce 
the influence of the experimental data in evaluating the dispersion integral. The amount of 
reduction increases with the order N of the polynomial Pn{s) that is being used to fit the 
weight function H(s) (N should not be chosen too large as will be explained later on). The 
second integral on the right hand side of Eq. (^) will be evaluated on a circle in the complex 
plane making use of identity (|5|) and state-of-the-art QCD input. In order to be explicit 
we evaluate the dispersion integral (0) between threshold and Si by the approximate sum 
rule 

12tt [ Sl lmIl(s)H(s)ds= f* R(s)(H(s) - P N (s))ds + 6m I Tl QCD (s)P N (s)ds. (6) 

Jso Jso J|s|=si 

It is generally believed that the results of perturbative QCD are valid on a circle of large 
radius in the complex plane, except possibly near the real axis. On the real axis instanton 
effects are likely to contribute significantly , invalidating the use of local duality, except 
for very large values of s. For example, the spectral function of the axial vector current 
extracted from r decay differs from the perturbative QCD prediction at s = up 
to a factor of three. In order to suppress the contributions to the contour integral close 
to the real axis we impose the condition that the polynomial P/v(s) vanishes for s = s\. 
Further conditions can be imposed on Pn{s) and are taylored to the specific problem at 
hand. For example, in the light quark case we impose the additional condition that the 
polynomial -Pjv(s) should coincide with the function H(s) in the p resonance region so as 



4 



to suppress the major contribution to the hadronic spectral function. In the heavy quark 
case to be discussed later on we will require that the difference (-Pjv(s) — H(s)) vanishes 
at the respective heavy quark thesholds where the quarkonium-resonances accumulate. In 
this way one can effectively suppress a large part of the experimental input needed for the 
evaluation of the dispersion integral (0). 

It may appear at first glance that, by increasing the degree N of the polynomial, 
the experimental input can be made arbitrary small. Although this would be true for 
the first integral on the right hand side of the sum rule Eq. (^J) one has to pay for this 
when evaluating the second integral where the QCD-input comes into play. The higher 
power terms in the polynomial put more and more emphasis on unknown higher order 
power suppressed terms, such as unknown higher order condensates for the light quarks or 
unaccounted for higher order terms in the mass expansion of the perturbative contributions 
of the heavy quarks. Moreover, since the coefficients of the polynomial approximations 
increase very rapidly with the degree N (with an alternating sign pattern for fixed N) the 
contributions of the higher order power suppressed terms become even more important as 
iV increases. The restriction on the degree of the polynomial approximation is correlated 
with the condition that the fitting range should not be too large. This is particularly 
true close to the poles of H(s) where H(s) shows rapid changes. For the light flavours we 
have therefore decided to move the point of coincidence from the p resonance to the point 
s = (1 GeV) 2 . This point is close enough to the p resonance to reduce the influence of the 
data substantially, and large enough to allow meaningful polynomial approximations to 
H(s) in this region of rapid change. The quality of the polynomial approximations in the 
light quark region can be judged by looking at Fig. 1(a) where polynomial approximations 
of different degrees N are plotted. We have checked on the consistency of our procedure by 
varying interval sizes and other constraints on the polynomial approximations and found 
no significant changes in the results. 

3 The light flavour contribution 

The first region corresponding to the light flavours u, d and s extends from 4m 2 to m% = 
(3.1 GeV) 2 , close to the charm threshold. In the region near the light quark threshold 
the function H(s) rises rapidly. This rise cannot be well reproduced by the polynomial 
approximation, but fortunately the annihilation cross section is very small near the light 
quark threshold. To obtain a reasonable low N polynomial approximation we therefore 
minimize (H(s) — Pn{s)) only in an interval from close to m 2 to m% apart from the 
additional conditions remarked on before. 

Turning to the QCD input for the light quarks we separately list the purely perturbative 
contributions and the condensate contributions. For massless quarks the imaginary part 
of the two-point function is known up to four loops in QCD perturbation theory ||. For 
the strange quark we include the 0(m 2 /s) power correction to three- loop order The 
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perturbative contribution to the current- current correlator reads ||, 
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(7) 



with k — 1, ki — 1.63982 and k 2 = 6.37101 (for the $ see below). We have denoted the 
unknown constant term in the four-loop contribution by C3. The constant non-logarithmic 
terms do not contribute to the circle integrals. The condensate contributions which we will 
refer to as the non-perturbative contributions are given by [TTD] 
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where we have included the m^/s 2 -contribution arising from the unit operator. In this 
expression we used the SU(3) colour factors Cf = 4/3, Ca = 3, Tp = 1/2 and L = 
ln(— fi 2 /s). The number of active flavours is denoted by n/. 

The result depends logarithmically on the renormalization scale \l and on the param- 
eters of the theory that are renormalized at the scale //. These are the strong coupling 
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constant, the quark masses and the condensates. As advocated in 12| , we implement the 
renormalization group improvement for the moments of the electromagnetic correlator by 
performing the integrations over the circle of radius s = s\ with constant parameters, 
i.e. they are renormalized at a fixed scale /i. Subsequently these parameters are evolved 
from this scale to fi 2 = si using the four-loop /3 function. In other words, we impose the 
renormalization group equation on the moments rather than on the correlator itself. This 
procedure is not only technically simpler but also avoids possible inconsistencies inherent 
to the usual approach where one applies the renormalization group to the correlator, ex- 
pands in powers of ln(s//i 2 ) and carries out the integration in the complex plane only at 
the end. In the present case the reference scale is given by A^g. 

For the coupling constant a s of the strong interaction we use the four-loop formula [13 
although a three-loop accuracy would be sufficient for the present application. We take 
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For the running quark mass we use the four-loop expression f[4" 
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and where 
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£(z) is Riemann's zeta function. Again we could have remained with three-loop accuracy 
in the running of the quark mass. 



4 The heavy flavour contributions 

Up to this point we have discussed in some detail how to deal with the contributions of 
the light quarks up to the charm quark threshold. Beyond charm quark threshold one has 
to incorporate the charm contribution in addition to the contribution of the light quark 
flavours. Further on when going beyond bottom quark threshold one has to include in 
addition the bottom contribution. Finally, above top quark threshold the top contribution 
has to be added. 

Let us begin by discussing the region between charm and bottom threshold s± — m% < 
s < s 2 . As concerns the polynomially weighted second integral in (0) the charm contribu- 
tion is obtained by the contour integration on the circle \s\ = s 2 whereas the light quark 
contribution is now obtained from the difference of the contour integrations at \s\ — s 2 and 
\s\ = Si, i.e. for the light quark contribution we now have 

f 2 -ImU(s)P N (s)ds = / U QCD {s)P N {s)ds U QCD {s)P N {s)ds. 

J Sl 7T 2m J\s\=s 2 ' 2%l J\s\=si 

(15) 

In the region between charm and bottom quark threshold the weight function H(s) is 
reasonably smooth and can be well approximated by polynomial functions of low degrees. 
The quality of the polynomial approximations is shown in Fig. 2(a) for iV = 1,4,5. As 
mentioned before we impose the condition Pn{s) = H{s) at s — si — m% in order to 
suppress the charmonium contribution with no restriction at the upper end s = s 2 = 
m\ = (9.46 GeV) 2 . In particular we do not set P/v(s) = at s = s 2 since instanton 
effects are expected to be negligible at these energies. One obtains a very good polynomial 
approximation H(s) starting with iV = 4 as Fig. 2(a) shows. 
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For the charm and bottom quark we use the QCD perturbative result in terms of an 
expansion in powers of m 2 / s to an even higher order than in the case of the strange quark. 
The current correlator is known up to three loops (i.e. 0(a s 2 )) up to order (m 2 /s) 6 [JO] 



For our purposes the 0((m 2 /s) 6 ) accuracy is quite sufficient since we need not go beyond 
N = 6. For example, for s = s 2 = m\ = (9.46 GeV) 2 the contribution of the 6t/i-order 
term in the charm mass expansion amounts to a tiny O(10 -8 ) effect. On the other hand 
the availability of the high power expansion of allows us to go to quite high degrees 
of the polynomial approximation. This reduces the influence of data on our results and 
further allows us to check on the consistency of our procedure by comparing results for 
different N. 

We would like to briefly comment on the mass dependence of the perturbative charm 
contribution. While the lowest order result for the spectral density decreases when the 
charm mass is increased the addition of the higher order terms leads to an increase of 
the spectral density with increasing mass [[□]]. Nevertheless the result of integrating the 
polynomially weighted current correlator IT p (s) over the circle decreases with mass. This 
can again be understood from analyticity: when doing the equivalent integration over the 
polynomially weighted imaginary part of IT p (s) the range of integration becomes smaller 
as the mass increases and the threshold moves up. 

The bottom quark contribution has its threshold at s 2 = m% = (9.46 GeV) 2 . There is 
no natural choice for the upper radius S3. On the one hand one would like to choose S3 
high enough so that one can safely start using local duality to evaluate the remaining part 
of the dispersion integral (§) from S3 to infinity by replacing R(s) by its QCD counter part. 
On the other hand the interval between S2 and S3 should not become too large to degrade 
the quality of the polynomial approximation to the weight function H(s). We strike a 
compromise between these two requirements and introduce two intervals of approximately 
equal spacing. The one interval extends from s 2 = m| = (9.46 GeV) 2 to s 3 = (30 GeV) 2 
and the second interval extends from S3 = (30 GeV) 2 to S4 = (40 GeV) 2 . The integrations 
in the two intervals are done as discussed before using contour integrals on circles or 
differences of these. In Fig 3(a) we show the quality of the polynomial fit for the interval 
s 2 < s < S3. A similar quality is obtained for the interval S3 < s < S4 but is not shown 
here. The remaining part of the dispersion integral (H) starting from S4 = (40 GeV) 2 up to 
infinity is done using local duality, i.e. we substitute perturbative QCD for R(s) in (fj). 



5 Results 

We begin the presentation of our results by discussing the light quark region in the interval 
from 4m 2 = (0.28 GeV) 2 to (3.1 GeV) 2 . At present the condensate power corrections are 
only known up to order 0(s~ 4 ) with any degree of confidence. Therefore the maximal 
degree of the polynomial approximation is restricted to lie below N = 3 in the light quark 
region. For polynomials of higher degree, by Cauchy's theorem, unknown higher dimension 
condensates would contribute to the circular integration. In principle one could use the 
present methods to obtain information on the higher dimension condensates. However, 
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degree 


without approximation 


iV= 1 


N = 2 


iV = 3 


N = A 


experiment |J 
perturbative 
nonperturbative 


75.688 


-19.396 
92.416 
-0.174 


16.268 
59.051 
-0.518 


11.370 
61.693 
-1.469 


10.313 
62.704 
-2.019 


total 


75.688 


72.846 


74.801 


71.594 


70.998 



Table 1: Contributions to Aa{j ad (Mz) in the light quark region for different degrees of the 
polynomial approximation. We detail the contributions coming from different parts of the 
sum rule Eq. (^). First row: experimental contribution (first term on the r.h.s. of (^|)). 
Second and third row: perturbative and nonperturbative condensate contributions (second 
term on the r.h.s. of ([])). 

this avenue will not be pursued in the present paper. 

In Fig. 1(b) we plot our results for the light quark contribution to Aa^ d (Mz) for the 
lowest energy interval ranging from 4m 2 = (0.28 GeV) 2 to (3.1 GeV) 2 . The horizontal line 
gives the experimental result according to the l.h.s. of the sum rule Eq. @ and the zig- 
zag lines give the results of the evaluation according to the r.h.s. of Eq. (|6|) for different 
polynomial degrees N (the normalization is that of Eq. @). The error margins in Fig. 1(b) 
result from theoretical errors in the evaluation of the r.h.s. of Eq. (|6]). They are given by 
the condensate errors, by the error on the strange quark mass and the QCD scale error 
to be detailed later on. The quality of the polynomial approximations in this interval is 
shown in Fig. 1(a). In Table [I] we list numerical values for the various contributions to 
Aa^ d (Mz) in this interval which appear in the evaluation of the l.h.s. and r.h.s. of the sum 
rule Eq. (^). We have separately listed the perturbative and nonperturbative contributions 
to the second term on the r.h.s. of Eq. We have included also some results for N > 3 
to indicate the tendency of the calculation when N becomes larger. 

Already for the linear approximation N=l the contribution of the data to the the 
r.h.s. of the sum rule Eq. (||) is reduced significantly. When the degree of the polynomial 
approximation is increased from one to three, the influence of the experimental input is 
reduced step by step. Unfortunately, the contribution of the poorly known condensates 
rises simultaneously. We take the standard values for the condensate terms and assign 
generous errors of 100% to them. We thus have 

(-GG) = (0.04 ± 0.04) GeV 4 , ct s (qq) 2 = (4 ± 4) • 10~ 4 GeV 6 . (16) 

7T 

For the errors coming from the uncertainty of the QCD scale we take 

= 380 ± 60 MeV (17) 

The errors resulting from the uncertainty in the QCD scale in different energy intervals 
are clearly correlated and will have to be added linearly in the end. We also include the 
error of the strange quark mass in the light quark region which is taken as 

m s (lGeV) = 200 ±60 MeV (18) 
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Table [I] shows that one obtains consistent results for different choices of polynomial 
degrees. The consistency of the different results is non-trivial since the final result arises 
from the sum of very different numbers for each polynomial degree. As a further check 
on the consistency of our approach we have subdivided the interval of our polynomial fit 
into two smaller intervals, namely from s = m 2 p to m^, where perturbative QCD is already 
applicable, and from s = vr? T to m\. The results are very similar to the former calculation 
but lead to slighty larger errors. 

Next we discuss the choice of the central value for the dispersive sum rule evaluation 
and the methodological error that we assign to it. The result of our sum rule evaluation 
in any given energy interval depends on the degree N of the polynomial approximation 
(see Figs. 1(b), 2(b) and 3(b)). We therefore choose a pair of neighbouring values of iV 
to determine our central value and its error. The central value is determined by the mean 
of the two sum rule values and the methodological error is given by the deviations from 
this mean. The choice of iV's is determined by the following two conflicting criteria. First 
N should be as big as possible so as to reduce the influence of the data. Second, there 
should be as little contribution from the poorly known condensates as possible. Latter 
criterion is only important for the first light quark energy interval. For the higher lying 
energy intervals we take median values of N. Our choices of neighbouring pairs of N in 
the various energy intervals are listed in Table |2|. 

As Figs. 2(a) and 3(a) show the polynomial approximations to the weight function 
become better and better as one is moving away from the lower pole of the weight function 
H(s) at s=0. This implies that the influence of the data on the sum rule evaluation 
becomes smaller and smaller as one is moving up in energy. This is quite apparent in 
Table 2 were we have listed the fractions of the experimental contribution to the sum 
rule for the different energy intervals. The data errors in the different energy intervals are 
multiplied by the same percentage figures and are added in quadrature to the final error. It 
is clear that the bulk of the experimental error comes from the lowest light quark interval 
while the contribution from the higher lying energy intervals are negligibly small. 

The contributions of the remaining energy intervals are collected in Table |2|. The large 
error for the second interval starting at charm threshold results mainly from the large error 
in the charm quark mass. For the charm and bottom quark masses we use the values 

fh c {m c ) = 1.4±0.2GeV, rh b (m b ) = 4.8 ± 0.3 GeV. (19) 

Summing up the contributions from the five flavours u, d, s, c and b our result for the 
hadronic contribution to the dispersion integral including the systematic error due to the 
dependence on A^jg (column 5 in Table 2) reads 

Aa^ d (M z ) = (277.6 ± 4.1) • 10" 4 . (20) 

In order to obtain the total result for a(Mz), we have to add the lepton and top contri- 
butions. Since we have nothing new to add to the calculation of these contributions we 
simply take their values from 0, who quote 

Aa^ ad (M z ) = (-0.70 ± 0.05) • 10" 4 , Aa lep (M z ) « 314.97 ■ 10~ 4 . (21) 
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interval 
for y/s 


values 
of iV 


data 
contribution 


contribution 

to Aa^ d (M z ) 


error 
due to Aj^g 


[0.28 GeV, 3.1 GeV] 
[3.1 GeV, 9.46 GeV] 
[9.46 GeV, 30 GeV] 
[30 GeV, 40 GeV] 
y/s > 40 GeV 


1,2 
3,4 
3,4 
3,4 


24% 
0.3% 
1.1% 
0.15% 


(73.9 ± 1.1) • 10" 4 
(69.5 ±3.0) • 10" 4 
(71.6 ±0.5) • 10" 4 
(19.93 ±0.01) • 10~ 4 
(42.67 ±0.09) • 10~ 4 


0.9 • 10" 4 
1.4- 10" 4 
0.06 • 10" 4 
0.02 • 10" 4 


total range 






(277.6 ±3.2) ■ 10" 4 


1.67- 10" 4 



Table 2: Contributions of different energy intervals to a\J d (Mz)- Second column: choice 
of neighbouring pairs of the polynomial degree N. Third column: fraction of the contribu- 
tion of experimental data [l|]. Fourth column: contribution to Aa^ d (M z ) with all errors 
included except for the systematic error due to the dependence on A^jg which is separately 
listed in the fifth column. 

Writing Aa(M z ) = Aa lep (M z ) + Aa had (M z ) our final result is (a(O)" 1 = 137.036) 

a(M z y l = a(0) _1 (l - Aa{M z )) = 128.925 ± 0.056. (22) 

6 Conclusion 

We have presented a new determination of the running fine structure constant at the scale 
of M z where we have made use of theoretical QCD results to reduce the contribution of 
experimental e + e~ annihilation data. Our calculations are in a sense complementary to 
the results of two recent papers H [5|]. These authors use QCD to replace data || or 
renormalize data || in regions where different data sets are mutually inconsistent with 
the aim of reducing the error estimate of the purely phenomenological calculations [|TJ . In 
this way inevitably some model dependence is introduced through the estimates of the 
systematic errors in the experimental data. In our analysis this model dependence arising 
from estimates of the systematic errors in the data is greatly reduced. In [|3| local duality is 
assumed to hold for the light flavours down to a scale of about m 2 with contrary evidence 
from r decay concerning this assumption. In [Bj data in the resonance region are rescaled 
by renormalizing the same set of data outside the resonance region to QCD, the assumption 
being that the normalization of the data did not change with energy and time over the 
long period of time while they were taken. In addition only three of the available eight 
data sets could be treated in this manner. In this way the authors of both papers succeed 
in substantially reducing the errors on a(M z ). 

In contrast to this, our philosophy is complementary. We want to arrive at a conser- 
vative result on a(M z ) and its error of which is as free of assumptions as possible. For 
instance, we use global duality at the large scale of Si = (3.1 GeV) 2 and, for all practical 
purposes, eliminate the effect of instantons by using duality only for combinations of poly- 
nomials which vanish at s\ = (3.1 GeV) 2 . We mention that our results are consistent with 
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those of H [5] . All three results are of great interest in the search for the Higgs boson. Our 
results should to be used if a conservative window for the Higgs mass is desired. 

We would like to close with the remark that all three recent calculations of at(Mz) 
should not deter experimentalists from remeasuring the e + e~ annihilation cross section 
more accurately in the low and intermediate energy region, as such data are absolutely 
essential for a precise value of a(M z ), unbiased by theory. 



Acknowledgements: We would like to thank F. Jegerlehner for providing us with the 
data set used in jjj]. We also want to thank R. Harlander for sending us the MATHEMAT- 
ICA output of the formulas presented in |TTJ. Further we would like to thank A.H. Hoang 



for correspondence and G. Quast for discussions on all experimental aspects of this work 
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Note added in proof: When preparing the final version of this paper for publication 
we received a preprint of M. Davier and A. Hocker on the same subject ( |hep-ph / 9805470]) 
using techniques similar to those described in this paper. 
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Figure Captions 



Fig. 1: (a) Weight function H(s) and polynomial approximations Pn{s) in the low- 
est energy interval 2m n < y/s < 3.1 GeV. The least square fit was done in 
the interval m p < y/s < 3.1 GeV with further constraints H(s) = Pn(s) at 
y/s = 1 GeV and Pjv(s) = at y/s = 3.1 GeV. The quality of the polynomial 
approximations are shown up to iV = 4. We use the scaled variable s/s\ 
for the polynomial approximation where s\ is the upper radius such that 
Pn( s / s i) is dimensionless. 

(b) Comparison of the l.h.s. and r.h.s. of the sum rule Eq. (Q) in the interval 

0. 28 GeV < y/s < 3.1 GeV. Dotted horizontal line: value of integrating the 

1. h.s. using experimental data including error bars fl[]. The points give the 
values of the r.h.s. integration for various orders N of the polynomial approx- 
imation. Straight line interpolations between the points are for illustration 
only. The dashed lines indicate the error estimate of our calculation. 

Fig. 2: (a) Weight function H(s) and polynomial approximations Pn{s) in the in- 
terval 3.1 GeV < y/s < 9.46 GeV with further constraint H(s) = -Pjv(s) at 
y/s = 3.1 GeV. Shown are the polynomial approximations for N = 1,4,5 
where we use the scaled variable s/s2 in the polynomial approximations, 
(b) Comparison of the l.h.s. and r.h.s. of the sum rule Eq. (|6|) in the interval 
3.1 GeV < y/s < 9.46 GeV. Dotted horizontal line: value of integrating the 
l.h.s. using experimental data including error bars [l|]. The points give the 
values of the r.h.s. integration for various orders N of the polynomial approx- 
imation. Straight line interpolations between the points are for illustration 
only. The dashed lines indicate the error estimate of our calculation. 

Fig. 3: (a) Weight function H(s) and polynomial approximations Pn{s) in the in- 
terval 9.46 GeV < y/s < 30 GeV with further constraint H(s) = Pn{s) at 
y/s = 9.46 GeV. Shown are the polynomial approximations for N — 1,4, 5 
where we use the scaled variable s/s2 in the polynomial approximation, 
(b) Comparison of the l.h.s. and r.h.s. of the sum rule Eq. @ in the interval 
9.46 GeV < y/s < 30 GeV. Dotted horizontal line: value of integrating the 
l.h.s. using experimental data including error bars [l|]. The points give the 
values of the r.h.s. integration for various orders N of the polynomial approx- 
imation. Straight line interpolations between the points are for illustration 
only. The dashed lines indicate the error estimate of our calculation. 



15 



20 
15 
10 
5 

5 



t^tt i i | iv i i i | i i i i | i i i i | i i i i | i i i i | i i i i | t i i r~i 

jn units 1 0~ 3 ) pole weight ; >- 

linear approx. 
quadratic appro/x- 
cubic approx. / - 
quartic approx'. - 




0.5 1 



.5 2 2.5 



Figure 1(a) 



5 3.5 4 
E (in GeV) 



90 



80 



i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i r 

theory (mean value) m s = 200 ± 60 M eV 

theory (error estimate) 

experimental mean value 

experimental error estimate 



'0 



fit from Vs = 1 GeV 
to Vs = 3.1 GeV 



60 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 







1 2 3 4 5 6 7 

degree of the polynomial approximation 



Figure 1(b) 



20 
15 
10 
5 

■5 



1 i i i | i i i rn i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i r 

pole weight 

linear approx. 

quartic approx. 

quintic approx. 

(in units 1 0" 3 ) 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 



1 2 3 4 5 6 



8 9 10 
E (in GeV) 



Figure 2(a) 



iO 



70 



60 



n i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i r 

theory (mean value) ryi c = 1 .4i0.2 GeV 

_ theory (error estimate) 

experimental mean value 

experimental error estimate 



fit from Vs = 3.1 GeV j 

t to Vs = 9.46 GeV : 
_i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 

1 2 3 4 5 6; 

degree of the polynomial approximation 



Figure 2(b) 



20 
15 
10 
5 

5 



n i i i | r^i r i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i r 

pole weight 

linear approx. 

quartic approx. 

quintic approx. 

(in units 1 0" 3 








J J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I Lj I I L_ 

5 10 15 20 25 30 35 40 

E (in GeV) 



Figure 3(a) 



90 



80 



i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i r 

theory (mean value) m b = 4.8 ± 0.3 GeV 

theory (error estimate) 

experimental mean value 

experimental error estimate 



'0 



fit from Vs = 9.46 GeV 
to Vs = 30 GeV 



60 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 







1 2 3 4 5 6 7 

degree of the polynomial approximation 



Figure 3(b) 



